Tight lower bound on geometric discord of bipartite states 
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We use singular value decomposition to derive a tight lower bound for geometric discord of arbitrary bipartite 
states. In a single shot this also leads to an upper bound of measurement induced non locality which in turn 
yields that for Werner and isotropic states the two measures coincide. We also emphasize that our lower bound 
is saturated for all 2 ® n states. Using this we show that both the generalized GHZ and W states of N qubits 
satisfy monogamy of geometric discord. Indeed, the same holds for all A'-qubit pure states which are equivalent 
to W states under stochastic local operations and classical communication. We show by giving an example that 
not all pure states of four or higher qubits satisfy monogamy. 
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Recent years have witnessed the emergence of some non 
classical correlations other than entanglement. Of them, the 
quantum discord is the most well studied and it indicates that 
separable states may possess quantumness which can be ex- 
ploited in various tasks e.g., state merging. There are different 
versions of quantum discord and their measures. However, al- 
most all measures are very difficult to calculate analytically, 
except the geometric discord (GD) introduced by Dakic et al. 
[1]. GD is defined as 



Dip) = mm\\p- X t 



(1) 



where Qo is the set of zero-discord states (i.e., classical- 
quantum states, given by PklftkXftkl ® Pk) an< i \\A\\ 2 = 
Tr(A 1 A) is the Frobenius or Hilbert-Schmidt norm. The au- 
thors in [1] have also calculated D for arbitrary 2-qubit states, 
using the explicit Bloch representation. This however poses a 
problem in generalizing the formula since the explicit Bloch 
representation is not known beyond 2-qubits (particularly, 
conditions for a vector v e R''" -1 to represent the Bloch vec- 
tor of a qu-rfit is not known for d > 3). So, this problem can 
not be solved analytically, in general. Fortunately, Luo and 
Fu have given an alternative description of GD in [2], via a 
minimization over all possible von Neumann measurements 
on p" 



D(p) = min ||p- 



n"(p)|| 2 



(2) 



and cast GD as the following optimization problem: 

Dip) = Tr(CC') - max Tr(ACC'A') (3) 

A 

where C = iCij) is an m 2 x n 2 matrix, given by the expansion 

p = Y l c i jX i ®Y j (4) 

in terms of orthonormal operators X/ e LiH a ), Yj e LiH h ) and 
A = (a*,) is an m x m 2 matrix given by 



aa = jmmxi = mxm 



(5) 
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for any orthonormal basis {\k}} of H a . Thus, the problem 
of determination of D reduces to finding the maximum of 
/(A) :=Tr(ACC'A') subject to the restriction in (5). Some 
effort has been directed towards this last part [3]. In this Brief 
Report, we derive a lower bound of GD for arbitrary states 
which will be shown to be saturated by all 2 ® « states. 

Another post-entanglement measure of quantum correla- 
tions is the measurement induced nonlocality (MIN), intro- 
duced by Luo and Fu [4]. The MIN is defined as somewhat 
dual to the GD, by 



Af(p) = max|b-n»|| 2 



(6) 



where the maximum is taken over the von Neumann measure- 
ments IT' = {Yl a k } which do not disturb p" locally, that is 



Y J ^p a K=p a 



(7) 



Thus, MIN is an indicator of the global effect on the whole 
system p ab caused by locally invariant measurement applied 
to one part, p a . In [4], the authors have calculated MIN for 
arbitrary pure states and 2 ® n mixed states. We will show 
that our lower bound of GD automatically reduces to the up- 
per bound of MIN derived therein. Using this bound, we show 
that the Werner and isotropic states have same amount of GD 
and MIN. These states are good candidates for maximally en- 
tangled states and have been studied frequently in literature. 

A tight lower bound on geometric Discord for arbitrary 
states: To solve the optimization problem in (3), it is helpful 
to fix the orthonormal bases {Xf], {Yj} and usually the follow- 
ing Bloch representation is considered: 

p = — \l m ® /„ + x'A <g> /„ + I m ® fX + V T u Ai ® Aj] (8) 
mn 1 1 

where A = (Ai, A2, ■ ■ ■ , A d i_ x )' with A, being the generators of 
SUid) for appropriate dimension d = m,n [5]. Comparing 
the two forms of p given by (4) and (8), we identify X\ 

J_/ 

VST"" 



Xit\ - ^j=Ai-uYj±i - -^Aj-i and 



C = 



y/mn 



s/mn 



(9) 
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Next we observe that the restriction (5) basically gives the fol- 
lowing three restrictions on A: 

e':M«ii)l''r(«ICi = "4=(1. L • • •. D (10a) 

ym 



in in in 



/,= ! 



(TrX,)r; 2 = (10b) 



the isometry condition AA' — I m (10c) 
and \k){k\ should be a legitimate pure state. (lOd) 

Before proceeding further, we note that the condition (lOd) 
means (a#)7Lj should be a coherence vector for all and as 
mentioned before, there is no known sufficient condition for 
it beyond R 3 . Thus, this constraint generically can not be im- 
plemented into the optimization problem for m > 3. So, for 
the the time being, let us ignore this constraint and optimize 
(maximize) f(A) with respect to the other constraints. Clearly, 
that would give us a lower bound of D(p). 

To incorporate (10a) into A, we write A — (e B) where B 
is any m x m 2 - 1 matrix subject to the restrictions (10b) and 
(10c). With these forms of A and C, we have 



f(A) = — 
mn 



l + 2 - 

n 



+ 2TiiB 



2 2 V2 
-x+— Ty 

n ym 



+Tr^fil-xx' + —TT'\B' 

\ m mn 



(11) 



Noting that xe' = -t=(x,x, . . . ,x), we have Tr(Bxe') = 

Y!k=i a k- x - 0, by (10b). Similarly, noting that Ty is a col- 
umn vector, we have Ti(BTye') = and hence the first trace 
term in (11) vanishes. So, we are left with only the second 
trace term. 

Writing A = (e B), we have from (10c), (e B)(e> B'J = I m , 
or ee' + BB' = /,„. Thus B must satisfy 



BB' = I m - ee' 



(12) 



This shows that the eigenvalues of BB' are 1 (with multiplicity 
m— 1) and (with e being an eigenvector). Let us choose an 
m x m 2 - 1 orthogonal matrix U having e as its last column. 
Then, every B satisfying (12) can be written as B = ULV, 
where V is an m 2 - 1 x m 2 - 1 orthogonal matrix and 2 is 
an m x m 2 - 1 diagonal matrix with diagonal (1, 1, . . . , 1, 0). 
Then defining G :- ( — xx' + —TT'\ for brevity, the last term 
in (1 1) becomes 



g(B) = Ti[BGB'] = Ti [UZV'GVL'U 1 ] 
= Tr [s'l/'t/EV'GV] = Tr [AV'GV] 



(13) 



where A := IfE = diag(/ m _i, m 2_ m ). This shows that max- 
imum of g(B) occurs when V'GV is a diagonal matrix whose 
diagonal entries are in non-increasing order. Since G is real 
symmetric, there always exists such an orthogonal V. Hence 
we have 



max 



(14) 



where X\ are the eigenvalues of G sorted in non-increasing 
order. Substituting this value of g(B) in (1 1), we get max /(A) 
which in turn gives the desired lower bound for GD from (3) 

as 



D(p)> — 
mn 



-IWi 2 + — ||r|| 2 -V^ 

m mn *■ — ' 



k=1 



(15) 



We note that this straightforward derivation uses singular 
value decomposition and does not require any upper bound 
for /(A). This is an important advantage because it directly 
shows what the minimum of g(B) should be (which would cor- 
responds to min/(A) and will be needed for deriving MIN). 
A lower bound of GD has been derived in [2] using only 
the isometry condition (10c). Since we have used more con- 
straints, undoubtedly our bound is sharper . 

Before applying this lower bound to solve some interest- 
ing related problems, let us show that this bound could be 
achieved by an infinite number of (collection of measurement- 
like) operators IF = {|&)(&|}, where each \k)(k\ is a hermitian, 
unit trace, but not necessarily positive operator. If all \k)(k\ sat- 
isfy (lOd), it would correspond to the optimal von Neumann 
measurement IF, which would yield the minimum of GD. We 
note that IT 1 = {\k){k\} where 

|*><*| = Y.auXi = -I m + -=8*1, *= 1,2,. ..m- 1 (16) 

v2 



and \m){m\ = I m -Y!k=\ \k){k\. Thus we need to determine only 
the first (m - 1) projections \k)(k\ and for this we should con- 
sider only the first m— 1 rows of B. So, denoting corresponding 
restrictions of B, e, U, E by B m _i,e m _i, U,„-i,A respectively, 
(12) reduces to B m ^B' m ^ = / m _j - (1 - l/m)e m _ie^_ r This 
in turn gives 



B m _i = (ai, a 2 , . . . , a m _i)' = C/ m _i AV 



(17) 



where U m -\ has e,„_i as its last column, 
A =diag(l, 1, . . ., 1, 1/ ym) and columns of V are the 
eigenvectors of G corresponding to eigenvalues A k . We note 
that different choice of U m -\ corresponds to different \k){k\ 
(though the set 11" may remain invariant). For a particular 
explicit representation, out of many choices for the rest of the 
columns, a particular one is to choose U m -\ as the Helmert 
matrix [6] which is given by (for clarity column vectors are 
not normalized) 



m-l 



1 1 


1 


1 .. 


1 1 A 


-1 


1 


1 .. 


1 1 





-2 


1 .. 


1 1 








-3 . 


1 1 








. 


1 1 


I 





. 


. — m + 2 1 , 



(18) 



Denoting the row vectors of U m -\ (with normalized columns) 
as r^, we have from (17), 



k=\ 



a k = r k V, k=\,2,...m-\ 



(19) 
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where = r£ ° (1, 1, . . . , 1, 1/ \Jm) 0°' is entrywise multipli- 
cation) and V is the m - 1 X m 2 - 1 left-upper block of V'. We 
emphasize that for m > 4, the choice U m -i is not unique e.g., 
for m = 5, £/,„_! can be taken as the standard 4x4 Hadamard 
matrix. 

Upper bound for MIN and its saturation by Werner and 
Isotropic states: To calculate MIN for a state p, we have to 
find minimum of Tr(ATT'A') where A has to satisfy an addi- 
tional constraint (7). As in the case of GD, ignoring (lOd) and 
(7) we would get an upper bound of MIN. Setting G - TT', 
we see that the required minimum is exactly the minimum of 
g(A) in (13). Hence just like (14), we have 



ming(A) = 2^1 



(20) 



k=\ 



Thus we have the following upper bound on MIN 



N{p)< 



. ,n-i . m -in 

mn K mrnr 

k~\ \ k-\ 



where a\ are the eigenvalues of TT' sorted in non- 
decreasing (non-increasing) order. We note that this upper 
bound is exactly the same as derived in [4]. If we set x = 0, 
the extra constraint (7) for MIN gets automatically satisfied. 
In addition, if all the eigenvalues are equal, the lower bound 
of Dip) in (15) and the upper bound of Nip) in (21) coincide. 
So, if one of the bounds saturates, necessarily we will have 
D = N. As an interesting consequence, we give the following 
two examples. The m x m dimensional Werner states 

m-z mz - 1 
P=— ; 1 + — F, ze[-l, 1] 



m — m m J — m 
withF := 2Hl*Xfl®l/><*l has 

(mz-l) 2 



D = N 



m(m — l)(m + l) 2 
For the m x m dimensional isotropic states 

1 - z m 2 z - 1 
P = -J— rl + —r-rmm, Z€[0,1] 
mr - 1 m - 1 

with \W) := 1/ 2r=i \ k ) ® \ k > we have 



D - N 



(m z z-iy 



m(m — V)(m + l) 2 



All 2®n states saturate our lower bound: Setting m — 2, we 
see from (18) the unique U\ is just 1 (seen as 1 x 1 matrix), and 
hence from (19), ai = 1/ V2vi. Then from (16), the unique 
measurement operators are given by 



|1)(1|= "(/2 + ViA) 



projectors and hence satisfy (lOd). Thus all 2 ®« states satu- 
rate our lower bound showing its tightness. We wish to men- 
tion that GD for these states have also been derived in [7], 
following the approach of [1]. 

One immediate consequence of the saturation of lower 
bound is that it readily gives GD for any qubit state. This 
in turn enables us to check monogamy relations etc. for qubit 
states. We will consider this case in the following paragraph. 

Geometric discord is monogamous for both generalized 
GHZ and W states of N qubits:- Recently many authors have 
studied monogamy property of different versions of quantum 
discord [8-10]. A correlation measure Q is said to be monoga- 
mous iff for any tripartite state p\ri (generalization to arbitrary 
state is straightforward) the following inequality holds 



Q{ Pll ) + Qipn) < <3(pi|2 3 ) 



(23) 



The authors of [8, 9] have shown that for (a specific measure 
of) quantum discord, all 3-qubit pure W-type states violate 
monogamy relation, while the GHZ-type states may or may 
not violate monogamy. Here we will show that the N-qubit 
generalized GHZ state \GGHZ) = a\00 ■ ■ -0> + Z>|11 • ■ • 1> and 
the generalized W states \GW) = £f =1 ^1001*0 ■ ■ • 0) both sat- 
isfy monogamy for GD 1 . 

Since GD is non-negative and any bipartite reduced density 
matrix (RDM) p\K of \GGHZ) is classical, Eq. (23) is au- 
tomatically satisfied for \GGHZ). Indeed, the relation holds 
for any arbitrary Schmidt-decomposable state 2 ~\fAi\ii ■ ■ ■ i). 
Thus, GD is monogamous for \GGHZ). 

In case of \GW), being pure, it should have a Schmidt de- 
composition over the cut 1|23 . . .JV and the Schmidt coeffi- 
cients (square-root of eigenvalues of p t ) are given by c\ and 

1 - c 2 Hence by the result of [4], the right hand side of Eq. 



4 



(23) becomes 2det(p0 = 2cf (1 - cj) = 2cf(c 2 + c 2 + ■ ■ ■ + c z N ). 
To evaluate the left hand side we note that the required bipar- 
tite RDMs are given by 



Pu 



1 





c 2 CiQ 

c\c k c\ 





(24) 



Expressing in Bloch form, we have x = (0, 0, 1 - 2c 2 ) and 
T =diag(2ciQ,2ciQ, 1 - 2c\ - 2c 1 ). Hence we have by our 
formula 

D(p lk ) = c 2 c 2 + imin{4c 2 c 2 ,(l-2c 2 ) 2 +(l-2c 2 -2c 2 ) 2 } (25) 



Thus, sum- 



Using minja, b] < a, this gives D(p\k) < 2cjCj 
ming over £'s our claim follows. 

One notable observation is that if we set all q,s equal 
(1/ yN), then Eq. (23) becomes an equality. This is quite 
remarkable, because it is known that the same relation holds 



|2)<2| = - {I 2 



MA) 



(22) 



Since vi (which is the eigenvector corresponding to the largest 
eigenvalue of G) has norm 1, both the operators in (22) are 



1 Just two days prior to this submission, in an interesting work Streltsov et.al. 
[13] have proven that all pure three qubit states satisfy monogamy of GD. 
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for the entanglement measure tangle t [ 1 1 ] , where the concept 
of monogamy appeared for the first time. 

We will now show that the result remains unchanged even 
if we add a term co|00 ■ ■ ■ 0} to \GW), i.e., if we consider class 
of states including all A^-qubit pure states which are equiva- 
lent to W states under stochastic local operations and classical 
communication (SLOCC) [12]. In this case the right hand side 
of Eq. (23) becomes 2det(pi) = 2c 2 (c 2 + c 2 + ■ ■ ■ + c 2 N ). To 
evaluate the left hand side, we note that each RDM p^ has 
x = (2cqCi, 0, 1 - 2c 2 ) and 



T = 



2ciQ 
2c\Ck 



2cqC\ 




2coQ 1 - 2c, 



?r 2 



Therefore eigenvalues of xx' + TT' are given by Ai = 4c 2 c 2 , 
A 2 ,3 =a± yfb where a = (1 - 2c]) 2 - 2c\{\ -c\-c\- c\) + 
4c 2 (c 2 + cf) and b = 8^[-(-l + 2c 2 + 2c]) 2 - 2(-l + 3c 2 + 
2c 2 l )c 2 k - 2c 4 k ) + a 2 . Noting that ||x|| 2 + ||r|| 2 = 8c 2 c 2 + 8c 2 c 2 + 
(1 - 2c]) 2 + 4c 2 c 2 + (1 - 2c\ - 2c 2 ) 2 := 8c 2 c 2 + c, we have 

||x|| 2 + ||n| 2 - maxU! , ^2,^3} < l|x|| 2 + ||r|| 2 - A 2 



= 8c 2 c 2 + c - (a + Vb) 



22 



< 8CTC7 



(26) 



where we have usedfr = (c-a) 2 + 32c 2 ) c 2 c^(l -Cg-c 2 -c 2 ) > 
(c - a) 2 . Hence D(p\k) < 2c 2 c\ and summing over k's the 
desired result follows. 

Due to this similarity with tangle it may be tempting to 
think that GD is also monogamous (at least) for all A^-qubit 



pure states. But GD, in contrast to tangle, is not monogamous 
for mixed states [13]. This indicates that may be GD is not 
monogamous for all pure states. To show this, let us consider 
the following A^-qubit pure state: 



m= Vpioo- ■■()>+ VWI + 1---D 



(27) 



where |+) = 1/V2(|0) + |1». For this state, we have 
D(pi\23...n) = 2det(pO = p(l - p), whereas D(p lk ) = 
l/2min[p 2 ,(l - p) 2 }. The state being symmetric in parties 
2, 3, . . . , N, monogamy relation (23) is satisfied iff 



N-\ 



min{p-,(l - pY) <p(l-p) 



(28) 



Clearly all 



N-l 



N+l N+l 



violate this relation. Thus not all pure states, beyond 3-qubits, 
satisfy monogamy of GD. 

To conclude, we have derived in a very simple way, a tight 
lower bound for geometric discord of arbitrary bipartite states 
which is saturated by all 2 ® n states. We have also shown 
that Werner and isotropic states have same amount of geomet- 
ric discord and measurement induced non locality. All pure 
A^-qubit generalized GHZ and W states are shown to satisfy 
monogamy of geometric discord. Giving an example we have 
shown that not all pure states of four or higher qubits satisfy 
monogamy of geometric discord. 

We would like to thank P. S. Joag for helpful discussions. 
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